Small RNA transcriptomic analysis during chronic viral hepatitis and cancer

Hepatitis B virus (HBV) and hepatitis C virus (HCV) are unrelated non-cytopathic viruses that infect the liver. Hundreds of millions of people worldwide are chronic carriers of HBV or HCV, and together these viral infections are responsible for ~80% of hepatocellular carcinoma (HCC). Numerous reports suggest that certain microRNAs (miRNAs, small non-coding RNAs that regulate gene expression post-transcriptionally) are important in the pathogenesis of viral hepatitis, although very few studies have assessed them comprehensively. Furthermore, no study of viral hepatitis or HCC has analyzed other classes of small RNAs, such as tRNA-derived RNAs (tDRs). tDRs have recently been suggested to play a role in human health, although they have never been previously assessed in primary human tissue. Our goal was to characterize thoroughly the small RNA transcriptome in chronic viral hepatitis and HCC. We used high-throughput sequencing to conduct an unbiased profiling of small (14-40 nts) RNAs in liver from Japanese subjects with advanced hepatitis B or C and paired HCC. We also analyzed gene expression data (microarray) and determined that pathways controlling cholesterol homeostasis were among the most significantly altered in chronic viral hepatitis and especially in HCC. An integrative analysis of miRNA and gene expression data predicted several candidate master miRNA regulators of each disease state, including miR-21 and miR-27. We validated in cell culture the roles of miR-21 and miR-27 as regulators of cholesterol synthesis, suggesting that these two miRNAs may be responsible in part for cholesterol imbalance in chronic viral hepatitis and HCC. We next sought to assess tDRs in chronic viral hepatitis and liver cancer. From the small RNA-seq, we observed that tDRs were abundant and significantly increased in chronic viral hepatitis. Remarkably, tDR abundance exceeded that of miRNAs in most infected non-cancerous tissue. In cancer tissue, the total abundance of tDRs was reduced and the relative abundance of individual tDRs was altered. Despite the potential relevance of tDRs to human health and disease, there is no standardized nomenclature, and no method has yet been developed to quantify tDRs from small RNA-seq datasets. To ameliorate these limitations, I developed a publically available method called tDRmapper.Doctor of Philosoph

tDRmapper: challenges and solutions to mapping, naming, and quantifying tRNA-derived RNAs from human small RNA-sequencing data

Abstract Background Small RNA-sequencing has revealed the diversity and high abundance of small RNAs derived from tRNAs, referred to as tRNA-derived RNAs. However, at present, there is no standardized nomenclature and there are no methods for accurate annotation and quantification of these small RNAs. tRNA-derived RNAs have unique features that limit the utility of conventional alignment tools and quantification methods. Results We describe here the challenges of mapping, naming, and quantifying tRNA-derived RNAs and present a novel method that addresses them, called tDRmapper. We then use tDRmapper to perform a comparative analysis of tRNA-derived RNA profiles across different human cell types and diseases. We found that (1) tRNA-derived RNA profiles can differ dramatically across different cell types and disease states, (2) that positions and types of chemical modifications of tRNA-derived RNAs vary by cell type and disease, and (3) that entirely different tRNA-derived RNA species can be produced from the same parental tRNA depending on the cell type. Conclusion tDRmappernot only provides a standardized nomenclature and quantification scheme, but also includes graphical visualization that facilitates the discovery of novel tRNA and tRNA-derived RNA biology

Prioritization of Genetic Variants in the microRNA Regulome as Functional Candidates in Genome-Wide Association Studies

Comprehensive analyses of results from genome-wide association studies (GWAS) have demonstrated that complex disease/trait-associated loci are enriched in gene regulatory regions of the genome. The search for causal regulatory variation has focused primarily on transcriptional elements, such as promoters and enhancers. microRNAs (miRNAs) are now widely appreciated as critical posttranscriptional regulators of gene expression and are thought to impart stability to biological systems. Naturally occurring genetic variation in the miRNA regulome is likely an important contributor to phenotypic variation in the human population. However, the extent to which polymorphic miRNA-mediated gene regulation underlies GWAS signals remains unclear. In this study, we have developed the most comprehensive bioinformatic analysis pipeline to date for cataloging and prioritizing variants in the miRNA regulome as functional candidates in GWAS. We highlight specific findings, including a variant in the promoter of the miRNA let-7 that may contribute to human height variation. We also provide a discussion of how our approach can be expanded in the future. Overall, we believe that the results of this study will be valuable for researchers interested in determining whether GWAS signals implicate the miRNA regulome in their disease/trait of interest

Differential hepatitis C virus RNA target site selection and host factor activities of naturally occurring miR-122 3′ variants

In addition to suppressing cellular gene expression, certain miRNAs potently facilitate replication of specific positive-strand RNA viruses. miR-122, a pro-viral hepatitis C virus (HCV) host factor, binds and recruits Ago2 to tandem sites (S1 and S2) near the 5΄ end of the HCV genome, stabilizing it and promoting its synthesis. HCV target site selection follows canonical miRNA rules, but how non-templated 3΄ miR-122 modifications impact this unconventional miRNA action is unknown. High-throughput sequencing revealed that a 22 nt miRNA with 3΄G (‘22–3΄G’) comprised <63% of total miR-122 in human liver, whereas other variants (23–3΄A, 23–3΄U, 21–3΄U) represented 11–17%. All loaded equivalently into Ago2, and when tested individually functioned comparably in suppressing gene expression. In contrast, 23–3΄A and 23–3΄U were more active than 22–3΄G in stabilizing HCV RNA and promoting its replication, whereas 21–3΄U was almost completely inactive. This lack of 21–3΄U HCV host factor activity correlated with reduced recruitment of Ago2 to the HCV S1 site. Additional experiments demonstrated strong preference for guanosine at nt 22 of miR-122. Our findings reveal the importance of non-templated 3΄ miR-122 modifications to its HCV host factor activity, and identify unexpected differences in miRNA requirements for host gene suppression versus RNA virus replication

Model of fibrolamellar hepatocellular carcinomas reveals striking enrichment in cancer stem cells

The aetiology of human fibrolamellar hepatocellular carcinomas (hFL-HCCs), cancers occurring increasingly in children to young adults, is poorly understood. We present a transplantable tumour line, maintained in immune-compromised mice, and validate it as a bona fide model of hFL-HCCs by multiple methods. RNA-seq analysis confirms the presence of a fusion transcript (DNAJB1-PRKACA) characteristic of hFL-HCC tumours. The hFL-HCC tumour line is highly enriched for cancer stem cells as indicated by limited dilution tumourigenicity assays, spheroid formation and flow cytometry. Immunohistochemistry on the hFL-HCC model, with parallel studies on 27 primary hFL-HCC tumours, provides robust evidence for expression of endodermal stem cell traits. Transcriptomic analyses of the tumour line and of multiple, normal hepatic lineage stages reveal a gene signature for hFL-HCCs closely resembling that of biliary tree stem cells-newly discovered precursors for liver and pancreas. This model offers unprecedented opportunities to investigate mechanisms underlying hFL-HCCs pathogenesis and potential therapies

A Gene Expression Panel is Accurate for Diagnosis and Monitoring Treatment of Eosinophilic Esophagitis in Adults

Eosinophilic esophagitis (EoE) can be difficult to diagnose. We aimed to evaluate whether a gene expression score could differentiate adult EoE cases from non-EoE controls and to determine whether scores normalized after treatment for EoE

microRNA-122 Abundance in Hepatocellular Carcinoma and Non-Tumor Liver Tissue from Japanese Patients with Persistent HCV versus HBV Infection

Mechanisms of hepatic carcinogenesis in chronic hepatitis B and hepatitis C are incompletely defined but often assumed to be similar and related to immune-mediated inflammation. Despite this, several studies hint at differences in expression of miR-122, a liver-specific microRNA with tumor suppressor properties, in hepatocellular carcinoma (HCC) associated with hepatitis B virus (HBV) versus hepatitis C virus (HCV) infection. Differences in the expression of miR-122 in these cancers would be of interest, as miR-122 is an essential host factor for HCV but not HBV replication. To determine whether the abundance of miR-122 in cancer tissue is influenced by the nature of the underlying virus infection, we measured miR-122 by qRT-PCR in paired tumor and non-tumor tissues from cohorts of HBV- and HCV-infected Japanese patients. miR-122 abundance was significantly reduced from normal in HBV-associated HCC, but not in liver cancer associated with HCV infection. This difference was independent of the degree of differentiation of the liver cancer. Surprisingly, we also found significant differences in miR-122 expression in non-tumor tissue, with miR-122 abundance reduced from normal in HCV- but not HBV-infected liver. Similar differences were observed in HCV- vs. HBV-infected chimpanzees. Among HCV-infected Japanese subjects, reductions in miR-122 abundance in non-tumor tissue were associated with a single nucleotide polymorphism near the IL28B gene that predicts poor response to interferon-based therapy (TG vs. TT genotype at rs8099917), and correlated negatively with the abundance of multiple interferon-stimulated gene transcripts. Reduced levels of miR-122 in chronic hepatitis C thus appear to be associated with endogenous interferon responses to the virus, while differences in miR-122 expression in HCV- versus HBV-associated HCC likely reflect virus-specific mechanisms contributing to carcinogenesis. The continued expression of miR-122 in HCV-associated HCC may signify an important role for HCV replication late in the progression to malignancy

Transcriptomic Analysis of Chronic Hepatitis B and C and Liver Cancer Reveals MicroRNA-Mediated Control of Cholesterol Synthesis Programs

ABSTRACT Chronic hepatitis B (CHB), chronic hepatitis C (CHC), and associated hepatocellular carcinoma (HCC) are characterized by cholesterol imbalance and dyslipidemia; however, the key regulatory drivers of these phenotypes are incompletely understood. Using gene expression microarrays and high-throughput sequencing of small RNAs, we performed integrative analysis of microRNA (miRNA) and gene expression in nonmalignant and matched cancer tissue samples from human subjects with CHB or CHC and HCC. We also carried out follow-up functional studies of specific miRNAs in a cell-based system. These studies led to four major findings. First, pathways affecting cholesterol homeostasis were among the most significantly overrepresented among genes dysregulated in chronic viral hepatitis and especially in tumor tissue. Second, for each disease state, specific miRNA signatures that included miRNAs not previously associated with chronic viral hepatitis, such as miR-1307 in CHC, were identified. Notably, a few miRNAs, including miR-27 and miR-224, were components of the miRNA signatures of all four disease states: CHB, CHC, CHB-associated HCC, and CHC-associated HCC. Third, using a statistical simulation method (miRHub) applied to the gene expression data, we identified candidate master miRNA regulators of pathways controlling cholesterol homeostasis in chronic viral hepatitis and HCC, including miR-21, miR-27, and miR-33. Last, we validated in human hepatoma cells that both miR-21 and miR-27 significantly repress cholesterol synthesis and that miR-27 does so in part through regulation of the gene that codes for the rate-limiting enzyme 3-hydroxy-3-methyl-glutaryl-coenzyme A (HMG-CoA) reductase ( HMGCR ). IMPORTANCE Hepatitis B virus (HBV) and hepatitis C virus (HCV) are phylogenetically unrelated hepatotropic viruses that persistently infect hundreds of millions of people world-wide, often leading to chronic liver disease and hepatocellular carcinoma (HCC). Chronic hepatitis B (CHB), chronic hepatitis C (CHC), and associated HCC often lead to cholesterol imbalance and dyslipidemia. However, the regulatory mechanisms underlying the dysregulation of lipid pathways in these disease states are incompletely understood. MicroRNAs (miRNAs) have emerged as critical modulators of lipid homeostasis. Here we use a blend of genomic, molecular, and biochemical strategies to identify key miRNAs that drive the lipid phenotypes of chronic viral hepatitis and HCC. These findings provide a panoramic view of the miRNA landscape in chronic viral hepatitis, which could contribute to the development of novel and more-effective miRNA-based therapeutic strategies

Small tRNA-derived RNAs are increased and more abundant than microRNAs in chronic hepatitis B and C

Persistent infections with hepatitis B virus (HBV) or hepatitis C virus (HCV) account for the majority of cases of hepatic cirrhosis and hepatocellular carcinoma (HCC) worldwide. Small, non-coding RNAs play important roles in virus-host interactions. We used high throughput sequencing to conduct an unbiased profiling of small (14-40 nts) RNAs in liver from Japanese subjects with advanced hepatitis B or C and hepatocellular carcinoma (HCC). Small RNAs derived from tRNAs, specifically 30–35 nucleotide-long 5′ tRNA-halves (5′ tRHs), were abundant in non-malignant liver and significantly increased in humans and chimpanzees with chronic viral hepatitis. 5′ tRH abundance exceeded microRNA abundance in most infected non-cancerous tissues. In contrast, in matched cancer tissue, 5′ tRH abundance was reduced, and relative abundance of individual 5′ tRHs was altered. In hepatitis B-associated HCC, 5′ tRH abundance correlated with expression of the tRNA-cleaving ribonuclease, angiogenin. These results demonstrate that tRHs are the most abundant small RNAs in chronically infected liver and that their abundance is altered in liver cancer

Model of fibrolamellar hepatocellular carcinomas reveals striking enrichment in cancer stem cells

The aetiology of human fibrolamellar hepatocellular carcinomas (hFL-HCCs), cancers occurring increasingly in children to young adults, is poorly understood. We present a transplantable tumour line, maintained in immune-compromised mice, and validate it as a bona fide model of hFL-HCCs by multiple methods. RNA-seq analysis confirms the presence of a fusion transcript (DNAJB1-PRKACA) characteristic of hFL-HCC tumours. The hFL-HCC tumour line is highly enriched for cancer stem cells as indicated by limited dilution tumourigenicity assays, spheroid formation and flow cytometry. Immunohistochemistry on the hFL-HCC model, with parallel studies on 27 primary hFL-HCC tumours, provides robust evidence for expression of endodermal stem cell traits. Transcriptomic analyses of the tumour line and of multiple, normal hepatic lineage stages reveal a gene signature for hFL-HCCs closely resembling that of biliary tree stem cells—newly discovered precursors for liver and pancreas. This model offers unprecedented opportunities to investigate mechanisms underlying hFL-HCCs pathogenesis and potential therapies